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1 . Introduction 

Observations of the £ Aurigae system early in the International 
Ultraviolet Explorer (IUE) program revealed an apparent UV excess 
shortward of !5G0A which was interpreted as the contribution of a hot dwarf 
companion (Hack and Selvelli 1979). Further observation, both before and 
during the ensuing eclipse cycle, revealed that the source of the UV excess 
was highly variable (Boehm, Ferluga and Hack 1984). However, as 
demonstrated by several researchers (Clarke 1981; Crivellari and Praderie 
1982; Basri, Clarke, and Haisch 1985, hereafter BCH) the IUE 
cross-disperser grating acts to redistribute a measurable percentage of 
the longer wavelength light into the range of the short wavelength prime 
(5WP) camera, causing errors in both line and continuum flux levels in the 
uncorrected spectra. 

There have been earlier attempts to correct the £ Aurigae IUE spectra 
for scattered light. Hack and Selvelli (1979) subtracted 500 flux units (FN) 
from the entire spectrum, treating the scattering as if it were a uniform 
background contamination, independent of wavelength. Parthasarathay and 
Lambert (1983) also recognized that scattered light might be important in 
establishing the properties, and indeed the very existence, of a hot 
companion. They used the descattering methods suggested by Clarke ( 1981 ) 
and Crivellari et al . ( 1980) and concluded that the signal shortward of 1250 
A was purely scattered light. They also suggested several alternatives to 
the hot companion model for the remaining UV excess shortward of 1600 A 
(see section 4) . 

It is clearly of crucial importance to our understanding of the nature 
of the UV excess that we carefully evaluate the level of scattered light 
contributing to it. In order to do so we have applied the BCH algorithm to a 
number of low dispersion IUE spectra of € Aurigae from very early 
pre-eclipse through the most recent post-eclipse epochs. As noted in BCH, 
scattered light is most significant when there is a great deal of contrast 
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between adjacent spectral regions. Indeed, € Aurigae has such a spectrum, 
being of the same spectral type as Canopus, the star chosen by 3CH as the 
best candidate for their scattered light study. 

We have found that the IUE spectra vary on timescales comparable to 
the optical photometry (Schmidtke 1985) in agreement with Ake (1985) and 
that they are indeed partially contaminated by scattered light. Even after 
correction for this instrumental effect, however, a significant, time 
dependent UV excess is still present. 


2. Data Reduction and the Descattering Procedure 

The BCH correction for scattered light in IUE spectra goes beyond the 
earlier work mentioned above in being more general and is therefore of 
wider applicability. The authors determined a “best fit" empirical scattering 
profile for the cross-disperser grating by convolving an assumed 
scattering profile with a “true’ Canopus spectrum and comparing the result 
to various IUE images of this FO lb-11 star. The “true" Canopus consisted of a 
Kurucz ( 1979) model (T?ff = 7500 K, log g = 2.0 and log A = 0.0) except in 
the spectral range 1420 < X < 3440, where TD-i and OAO II fluxes were 
used. The adopted scattering profile was derived from results of tests on a 
replica of a grating made from the same master as that used to produce the 
one actually flown aboard IUE (Mount and Fastie 1978). The profile was 
modified to include scattering from the cores of emission lines and extended 
by adding exponentially decreasing wings. Four free parameters are needed 
to describe these modifications to the Mount and Fastie profile. The 
interested reader is referred to the SCH paper for more details of these 
parameters and the errors associated with their determination. 

Once the best fit profile is determined it can be deconvolved from any 
IUE spectrum to return a “descattered" spectrum, and this is the approach 
that we used. The default, normalized scattering profile assumes a level 
of 4% scattering from a bin at the reference wavelength (arbitrarily 
selected at the midpoint of the combined short and long wavelength range) 
as one of the four parameters. That light is then redistributed to all other 
wavelengths in accordance with the shape of the adopted profile. For other 
than the reference bin the efficiency scales as (X r of/X) 2 . Hence, a value of 
4% at 2450 implies a 16% scattering efficiency in tne region of the Lyman a 
line. 


Because the light scattered into the range of the 3WP camera originates 
from longer wavelengths, we have merged all SWP spectra with LWR or LWP 
observations taken on the same day. The spectra input into our subsequent 
procedures are absolutely calibrated fluxes from the SWP camera for X s. 
1950 A and from the LWR or LWP cameras at X i 1950 A, yielding an 
effective wavelength coverage of 1150 i X i 3350 A. We have applied the 
absolute calibration of Bohlin and Holm (1980) to the SWP and LWR spectra 
and that of Cassatella and Harris (1982) to the LWP images. In addition, we 
have used the correction suggested by Holm, et al. (1982) for errors due 
to variations in the camera head amplifier temperature. 
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The BCH procedure was originally designed to work with input data in 5 
A steps, wider bins being inappropriate for spectral fine studies. Our 
primary concern in this work is the short wavelength continuum ; hence we 
are not similarly restricted. Moreover, the signal/noise is quite poor in 
certain spectral intervals of key importance, i. e. , 1150 x X x 1550 A (the 
region which purportedly shows a UV excess but which has the smallest 
intrinsic SWP flux) and 2000 x Xx 2400 ( the 'nearest* region of the long 
wavelength spectrum and hence the most sensitive to the actual shape of the 
redistribution function for scattering into the short wavelength region). 
Therefore, we averaged the data in 50 A wide bins, appropriately weighted 
to eliminate reseaux and other flagged points of questionable quality, and 
modified the original BCH code so as to allow input data of any reasonable 
bin size. 

These merged, absolutely calibrated, binned spectra were then passed 
to the descattering algorithm. In all cases we have used the 'best fit' 
parameters determined by BCH for the scattering profile, motivated in part 
by the similarity in spectral type of Canopus and e Aurigae. We have 
investigated the effect of a twofold reduction in the assumed scattering 
efficiency, however, in light of some evidence that the 4% value might be too 
large (c.f. , see discussion below). 


3. Results and Discussion 

Figure 1 depicts the wavelength and time dependence of the scattered 
light found by the 3CH procedure for spectra obtained throughout the 
eclipse. The quantity plotted is the scattered light fraction, defined as 

R(A,$) = 100 X (F°(M>) - F c (A,<£>))/F°(A,<f>) 

where F°(A,$) = the flux at wavelength A and eclipse phase <£, for the 

uncorrected (input) spectrum, and 

F C (A,$) = the flux at wavelength A and eclipse phase for the 
corrected (output) spectrum. 

A quadratic interpolation was applied to the data along the phase axis in 
order to display non-uniformly spaced points on a uniform scale. In most 
cases the scattered light fraction among the longer wavelength bins is 
insignificant, whereas for A x 1550 A R(A,<£) increases rapidly. One 
exception is the region AA 2000-2400 which, as mentioned above, is the 
least sensitive portion of the long wavelength cameras. 

Without even considering the vertical axis of Figure 1 we can easily 
detect an interesting facet of the phase dependence of R(A,$) in the short 
wavelength bins. Not surprisingly, the figure implies that the amount of 
scattered light is greatest when the star is brightest, i.e., before first 
contact, after fourth contact and during the mid-totality brightness 
enhancement. Since we would not expect an uneclipsed hot secondary to 
show the same time dependence as the eclipsed primary we interpret this 
result, as an independent verification that we are indeed dealing with 
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scattered light. Now if we consider the values along the vertical axis we 
see that up to 20% of the flux at the short wavelengths is scattered light. 
Two important questions then arise: 1) how reliable are the assumptions 
that set the scale of the vertical axis? and 2) what about the flux that is not 
scattered light? We shall consider each question in turn. 

In his application of the BCH algorithm to the spectra of comets 
Feldman (private communication) derived a scattering profile from the 
shape of an overexposed Lyman a emisssion line. He found an efficiency of 
scattering from line center of about 2%, compared to the 4% found by 3CH. 
This was partially compensated by a larger fractional scattering into the 
near wings, important for emission line studies but not for the continuum 
case under study here. In Figure 2 we show the results of a test 
comparison of the two cases, 4% and 2% scattering efficiency, 'for our most 
recent post-eclipse observation. Although the fractional scattered light is 
also reduced by the expected factor of two, we note that there is still 
agreement in the regions of overlapping uncertainty (± 1 a). Hence, it 
seems reasonable to consider the 2% case as a lower limit to the actual 
values of R(A,$), for this and all other epochs. 

A rough method for estimating the amount of scattered light in an IUE 
spectrum has recently beer, suggested by Imhoff ( 1 985 ) . Her empirically 
derived procedure determines the number of DN per second scattered into 
the 5WP camera range from the long wavelength spectrum at 2400 A. Since 
the BCH procedure sums the contribution from the entire long wavelength 
spectrum it might be expected that the values of R(A,$) derived from it 
would be larger than those determined from Imhoff’s 'rule of thumb'. 
However, it is evident from Figure 2 that this is not the case. Even with the 
rather large uncertainties Imhoffs approximation predicts a larger 
fractional scattered light than that found using the BCH algorithm. This 
suggests either that the 'rule of thumb* procedure overestimates the 
scattered light or that the profile adopted by BCH leads to a considerable 
underestimate. The latter seems unlikely, since the profile was determined 
from a best fit to a star of the same spectral type as £ Aurigae. Lacking 
further information as to the validity of the assumptions that were made in 
constructing the BCH profile (i.e., do the far wings really fall off 
exponentially?) we shall henceforth interpret the larger values derived 
from the Imhoff procedure as convenient upper limits to the full BCH 
treatment. We shall use the upper limit results in the later figures, in lieu 
of error bars, to represent the maximum correction for scattered light. 

Let us return to the second question posed in the discussion of Figure 
1. Having subtracted away the scattered light contribution in each 
wavelength bin for each of the spectra, we are nowin a position to examine 
the intrinsic UV flux r and note how it behaves as a function of 'eclipse 
phase*. Units of eclipse phase are defined such that first contact occurs at 
<fr= Q.QQ and last contact at 1.00, based on dates predicted by Gyldenkerne 
(1970) for the eclipse in optical wavelengths. In these units second 
contact occurs at = 0.200 and. third contact at 0.790. 

In Figure 3 we present eclipse light curves, both before and after the 
scattered light correction, for three representative wavelength bins, 1350, 
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1450, and 1600 A. In all cases we plot magnitude versus $, defined as 

m N ($) =-2.5 log F x (<fr) -21.1 

where the constant 21.1 refers to the absolute energy calibration of Vega 
derived by Hayes and Latham (1975). One can easily show that the 
difference between the light curve corrected for scattered light and the 
uncorrected curve is related to the R values discussed above by 

Sm^($) = m^ c (<£) - m^°(<l>) = -2.5 log- ( I - R(X,<£) ) 
where R is used here as a fraction, rather than a percentage. 

Inspection of Figure s demonstrates a number of interesting results. 
First, it can be easily seen that, with increasing wavelength, the corrected 
light curves fall closer and closer to the uncorrected curves. For X > 1700A 
the curves are virtually indistinguishable. This trend is entirely consistent 
with the result, noted in Figure 1, that the percent scattered light is 
greatest for the short wavelengths in the SWP camera and approaches zero 
at the longer wavelengths. Second, we note that the corrected and 
uncorrected light curves differ most at the "brightest" phases, i . e. , 
before and after the eclipse itself and during those phases associated with 
the central brightening feature (for those wavelengths where there is still a 
difference), again in accord with the results of Figure 1. We see that the 
effect of scattered light is to 'fill in' the curve, making it shallower than 
it should be at the shortest wavelengths. Most notably, it is quite 
disputable whether or not the curve at 1350 A is an eclipse curve, 
especially considering the large uncertainties concomitant with the low 
intrinsic flux levels. 

To address directly the question of the reality of a UV excess in the £ 
Aurigae system we examine the spectra after we have removed the 
scattered light component. Does evidence for an apparent hot component 
remain? The answer seems to be "yes’. For several of our observations a 
modest percentage of the UV excess is accounted for once we subtract away 
the scattered light. In others, almost none of the excess is removed. 
Figure 4 shows a typical example of the former case, in which subtracting 
even the upper limit of scattered light does not significantly reduce the UV 
excess at wavelengths below 1550 A. The format of Figure 4 is deliberately 
similar to that used by Hack and Selvelli (.1979), the first report of an 
excess in IUE spectra of € Aurigae to appear in the literature. 

When spectra (corrected for scattered light) from several different 
dates throughout the eclipse are plotted in a manner similar to that used in 
Figure 4, large variations are seen in slope of the continuum shortward of 
1600 A. Figure 5 is a schematic illustration of the short wavelength 
variations observed. The fluxes in Figure 5 have all been arbitrarily 
adjusted to yield the same value at 3300 A as the model atmosphere of 
Kurucz ( 1979) used in Figure 4, a procedure which presumably removes the 
eclipse time dependence from the data. The time variations that remain are 
associated either with the eclipse-independent fluctuations of the hot 
component itself, with gaps or tunnels in the occulting secondary (Wilson 
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1971; Aitner et al. 1984), or with the cepheid-like pulsations of the 
primary, which can be quite large at short wavelengths (Ake and Simon 
1984). Figure 6 shows the data at 1300, 1400 and 1500 A, adjusted in the 
same manner, plotted versus •time. There is some suggestion of cyclic 
behavior in the curves, and a detailed study of this time dependence is 
presently underway. 

In early February of 1984 the ultraviolet spectrometer (UV5) aboard 
the Voyager I spacecraft observed £ Aurigae (R. Polidan, private 
communication). The spectrum obtained is consistent with a star of 
spectral type 85; however, since a known B5 star is .clo.se enough to £ 
Aurigae to have been included in the 0. 1 degree by 0.9 degree aperture, the 
identification is suspect. An additional UVS observation is planned, from 
which it is hoped that the uncertainty in the identification of the UV source 
will be resolved. Ake and Simon (1984) noted that the line spectrum of £ 
Aurigae shortward of 1400 A does not match that of a 85 star. They 
concluded from this that whatever the hot component is, it is not a hot 
star, which suggests that the flux observed by the UVS instrument was due 
to the known B5 star. Nevertheless, the UVS data seems not to be 
inconsistent with the fluxes in the short wavelength !UE spectrum obtained in 
mid-February of' i984. We will follow with great interest future UVS 
observations of £ Aurigae, with the intent being -to compare the results with 
future 1UE data andU, B, and V photometry. 


4. Summary 

As a result of this work we have found that light scattered from the 
longer wavelengths constitutes a small but non-negligible, wavelength and 
time dependent fraction of the measured flux in the far UV. We have not 
been able to unambiguously rule out the reality of the UV excess. However, 
we note that there are still uncertainties in the assumed scattering profile. 
New measurements of the scattering properties of the cross-disperser 
grating are planned in order to verify the results of Mount and Fastie ( 1978) 
and extend the wavelength coverage into the far wings of the profile (8. 
Woodgate, private communication). The results of these measurements will 
no doubt reduce some of these uncertainties. For the present, we feel that 
the 3CH approach is a significant improvement over the methods heretofore 
available for the treatment of scattered light in IUE spectra. 


1UE imposed constraints on the particle size and temperature of the 
disk surrounding the. secondary were first discussed by Castel 11 et a!. 
(1982) and Chapman et al. (1983). Later infrared results confirmed the 
prediction that the disk temperature had to be significantly lower than that 
of the supergiant (Backman et al. 1984). The apparent contradiction 
inherent in the coexistence of an ultraviolet excess and an infrared 
temperature around 500 K was one of the primary motivations for 
undertaking this study. Future observations, such as those possible with 
the Voyager UVS and other instruments, will undoubtedly gain us new insight 
into the problem. Of particular interest, in the long term, is observation of 
the € Aurigae secondary at quadrature and, later still, in eclipse. 
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Meanwhile, we look to theoretical models which allow both a UV excess and 
a cold disk, such as star spots on- the surface of the primary 
(Parthasarathy and Lambert 1983), or polar flows associated with a double 
star embedded in the disk (Lissauer and Backman 1985). 

We acknowledge the assistance of the 1UE staff in acquiring the spectra 
of € Aurigae throughout the many months of the eclipse cycle, and of the 
staff of the Regional Data Analysis Facility at Goddard Space Flight Center, 
whose helpful advice in reducing the data was an invaluable resource. We 
are grateful to Dr. G. Basri for providing the FORTRAN code of the BCH 
procedure and for taking the time to clarify some of the finer points 
involved in its application. 
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Figure I 

The relative percentage of scattered light in each wavelength bin 
for each observation (data is interpolated to give uniform spacing 
m the time coordinate) The tic marks along the phase axis refer 
to the contact points of the eclipse Rn apparent trend toward 
larger amounts of scattered light during the uneclipsed <i e , 
brighter) phases can be seen, as wet! as during the m.d-eclipse 
brightness enhancement apparent in the light curves This behavior 
at the shortest woveiengths is expected if the signal is scattered 


Rule of thumb method 

8CH mi th 4% efficiency 
at A re f 

BCH with 2 % ef f i c i ency 
at A 

Error bars represent t O 
of uncertainty 
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figure 2 

Relative percent scattered light as a function of wavelength for 
a typical IUE observation of Epsilon Aurigae < <t> e « t 434 ) The 
effect of varying the scattering efficiency parameter is shown 
Although the two cases give percentages which differ by a factor 
Of two, thm overlappiryj uncer taint i es show that this pora*«Ur if 
not critical, at least for the case of continue** scattering Rt 
t ts upper limit the 4* cose approaches the rough estimate d«rived 
using Imhoff'S (1983) “rule of thumb' 
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Eclipse Phase 

Figure 3 

Eclipse light curves for three different wavelength pi ns In 
each case the filled dots represent the uncorrected data while 
the plus signs or# the som data corrected for scattered I ight 
using the BCH program with the default parameters, and the 
dashed I me is the curve corrected using the “rule of thumb' 
method fit wove lengths longward of 1700 fi the three curves ore 
virtually indist inguishabie 



fi typical example of the partial reduction m the apparent UU 
excess ofter subtracting the scattered tight component The 
thin line IS a Kurucz c ig?Q > model stellar atmosphere with 
Teff • 7500 K , log g - 1 0, and log fl * 0 0, adjusted to th« 
distance and angular diameter of Canopus The filled squares 
are the TO- 1 and 0fi0-2 data for Canopus, scaled as in Hack i 
Set veil i <1Q?9> The spectrum of Epsilon fiurigae well outside 
of eclipse is shown by the remaining curves O pen circles no 
correction for scattered light, heav y I > ne correction using 
the BCH code with default paramters, do t ted l i ne correction 
using lnihoff s rule of thump estimate 
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Figure 5 

fl demons tro t i on of the ti»e variation of the short wovelenqth 
spectra with variations due to the eclipse “filtered out" by 
adjusting all the continue to the value of the Kuruez (1979) 
model (open circles; same model as used in the previous figure) 
at 3300 A ftt the extreme short move length end, the variations 
range over almost 2 magnitudes 


Julian Date - 2,440,000 

5200 5400 5600 5800 6000 6200 



Logarithm of the Epsilon Rur i gae fluxes, corrected for scattered 
light and adjusted to the same value at 3300 A, as a function of 
time By scaling the fluxes m the near UU the time variation due 
to the eclipse is suppressed, leaving evidence of another time 
dependence, associated etther with the UU excess itself or with 
the cepheid-like pulsations of the primary Open circles 1300 ft, 
as ter i sks 1400 ft, filled squares 1500 ft 
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